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ABSTRACT 

Major mergers of disk galaxies are thought to be a substantial driver in galaxy evo- 
lution. To trace the fraction and the rate galaxies are in mergers over cosmic times, 
several observational techniques, including morphological selection criteria, have been 
developed over the last decade. We apply this morphological selection of mergers to 
21 cm radio emission line (Hi) column density images of spiral galaxies in nearby sur- 
veys. In this paper, we investigate how long a 1:1 merger is visible in Hi from N-body 
simulations. 

We evaluate the merger visibility times for selection criteria based on four pa- 
rameters: Concentration, Asymmetry, M20 and the Gini parameter of second order 
moment of the flux distribution (Gm)- Of three selection criteria used in the litera- 
ture, one based on Concentration and M20 works well for the H I perspective with a 
merger time scale of 0.4 Gyr. Of the three selection criteria defined in our previous 
paper, the Gm performs well and cleanly selects mergers for 0.69 Gyr. The other two 
criteria (A-M20 and C-M20), select isolated disks as well, but perform best for face-on, 
gas-rich disks (T mgr ~ 1 Gyr). The different visibility scales can be combined with 
the selected fractions of galaxies in any large H 1 survey to obtain merger rates in the 
nearby Universe. All-sky surveys such as WALLABY with ASKAP and the Medium 
Deep Survey with the APETIF instrument on Westerbork are set to revolutionize our 
^ I perspective on neutral hydrogen and will provide an accurate measure of the merger 

fraction and rate of the present epoch. 
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1 INTRODUCTION 

The merger rate of galaxies over cosmic times is one of 
the big outstanding questions in the evolution of galax- 
ies. Two recent major observational efforts seek to address 
the r ate of galaxy mergers, one using close galaxy pairs 
(e.g., iPatton et ai1ll997l : jLe Fevre et alJl200Ch. and one us- 
ing q uantified morphology (e.g.. IConselicell2003l ; lLotz et all 
l2004h . Both give a fraction of the population at a given red- 
shift that is merging (/ m ) which then can be converted to 
a merger rate as soon as one knows how long a merger is 
identifiable as either a galaxy pair or morphologically dis- 
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turbed galaxy. Hence, both methods need a timescale for 
which mergers are identified as such. In this series of pa- 
pers, we focus on the morphological approach, specifically 
the morphology of galaxies in the 21 cm emission line (Hi). 
Despite the generally lower resolution of the H I maps com- 
pared to optical images, we argue that the morphological 
signature of a merger event can be equally or better observed 
because (a) the H I disk extends well beyond the stellar disk 
(i.e., offers the same number of resolution elements as the 
optical disk) , and (b) the atomic gas is disturbed well before 
the stellar disk. Therefore, Hi morphology looks promising 
as an alternate tracer of mergers, provided representative 
enough volumes are surveyed. 

A popular measure of the merger signal in galaxy mor- 
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phology uses the Concentration - Asymmetry-Smoothness 
parameters (CAS, IConselicel |2003[). but addition al param- 
eters such as Gini (G) and M20 ( Lotz et al.ll2004l ) are used 
extensively as well. The merger signal is usually found in the 
restframe ultraviolet or B-band, calibrated by galaxies in the 
local Universe. The timescale on which a merger is identi- 
fiable is found from CAS or G and M20 measurements of 
N-body sim ulations of mergers l|Conselicej|2006l ; lLotz et al.l 
2008L l2010al lbh or alternatively from the decline in identified 
mergers between two nearby re dshift bins, assu ming only 
passive evolution in merger rate (|Conselicd 12009). Both ap- 
proaches give consistent estimates of the merger time-scale 
of about a Gyr or less, allowing for estimates of the typical 
number of mergers a massive galaxy undergoes during its 
lifetime (|Conselicel[200l . 120091 : lLotz et al.ll2008l ). A limita- 
tion is that morphological ide ntification of merg ers is most 
effective for gas-rich galaxies l|Lotz et al. I l2010al ) and varies 
with t he mass ratio of the merging galaxies (e.g-. lLotz et alJ 
l2010bh . 

All previous work on quantified galaxy morphology as 
a merger tracer ha focused on rest-frame ultraviolet or the 
blue side of optical. The benefits are clear for this technique 
as the higher redshift galaxies can be observed with the 
Hubble Space Telescope at these wavelengths at comparable 
spatial resolution to the local reference samples observed by 
GALEX or the Sloan Digital Sky Survey (SDSS). Star for- 
mation triggered by the merger increases surface brightness, 
easing their identification at higher redshift. 

However, in the coming decade, a new window on 
gas-rich mergers will be opening up with the commission- 
ing of new radio telescopes and instrument s: the South 
African Karoo Array Telescope (Me erKAT; Ijonasl 120071 : 
iBooth et aill2009l ; Ide Blok et aj1|2009T), the Australian SKA 
Pathfi nd er (A SKAP; I Johnstonl 120071 Ijohnston et~all l200l 
2008allbl. 120091 ) . the Extended Very Large Array (EVLA; 
Napier 2006) and the APERture Tile In Focus instrument 



(APERTIF: IVerheiien et al.ll2008l ; lOosterloo et al.ll2009l ) on 

the Westerbork Synthesis Radio Telescope (WSRT). Ulti- 
mately, this investment in radio observat ories will culminate 
in th e Square Kilometer Array (SKA; ICarilli fc Rawlinga 
120041 ). These observatories will extend the 21 cm line ob- 
servations of galaxies to high redshift and low column den- 
sities. Specifically, the effective all-sky survey of the nearby 
Universe with ASKAP (the WALLABY project, Koribalski 
et al. in preparation) and the Medium Deep Survey with 
WSRT/APERTIF spans an ideal volume for a local merger 
fraction and rate measurement. The Hi disk's morphology 
can be an alternative tracer of the merger fraction and con- 
sequently merger rate. There is already ample anecdotal 
evidence for the sensitivity of Hi morphol ogy to interac- 
tion (e.g., "The Hi Rogue Galaxy Gallery" . iHibbard et al.l 
120011 . http://www . nrao . edu/ astrores/HIrogues/ ) , and 
more quantified relations betw een de nsity and H I morphol- 
ogy exist (e.g.. iBouchard et al.ll2009T ). 

In this series of papers, we explore how well the param- 
eterized Hi morphology trace mergers in two local Hi sur- 
veys: t he H 1 Nearby Galaxy Survey ( THINGS JWalter et all 
200881) and the Westerb o rk Hi Spiral Projec t (WHISP, 



van der Hulst et al l l200ll : Ivan der Hulstl |2002| ). We have 
found that the H I column density maps - while at typically 
lower spatial resolution than the optical wavelengths- are 
just as sensitive if not more so to the effects of a merger, es- 



pecially when expressed in the morphologic al qualifiers cus- 
tomarily use d in optical or UV classification (|Holwerda et al.l 
I2009l . l2011b1 ). For the first time, it may be possible to select 
ongoing merging galaxies from their Hi morphology with- 
out any need for a human observer. This will greatly simplify 
the classification of galaxies and mergers in the upcoming 
all-sky H I surveys. 

In a subset of the WHISP sample we showed that there 
is a simple cut in morphological parameter space that di- 
vides th e merging disks from t he general non-merging pop- 
ulation (|Holwerda et al.ll201ld ). The question remains, how- 
ever, how long a merger is identifiable in this new window be- 
fore we can convert an observed me rger fractions in WHIS P 
to a merger rate (our fifth paper, iHolwerda et all l2011al ). 
In this paper, we explore the merger visibility timescale 
using N-body simulations of L* disks with cold gas disks, 
both evolving passively and undergoing a violent 1:1 merger 
with a similar disk. The organisation of this paper is as fol- 
lows: section [2] discusses the morphological parameters we 
use briefly, section [3] describes the simulations we used, how 
the H 1 map was generated, and the limitations of these sim- 
ulations and the Hi perspective. Section[4]is a discussion of 
our results from the N-body simulation, and section [5] our 
conclusions. 



2 MORPHOLOGY 

In this series of papers we use the CAS sys- 
tem (Conce n tratio n- Asymmetry-Sm o othnes s) from 



iBershadv et~ai1 (|2000l ): IConselice et~ai1 (|2000[) ; IConselicd 
( 2003t ). the Gini /Man system from lLotz et all (|2004l ) and 
our own parame ter Gm, the Gini para meter of the second 
order moment (|Holwerda et al.l l201ld ). Concentration is 
defined as: 



C = 5 



/rgo 
Vr 2 n 



(1) 



where r% is the radius which includes that percentage of 
the intensity of the object. In an image with n pixels with 
intensities at pixel position (i, j) and hso{i, j) is the 

value of the pixel in the rotated image, Asymmetry is: 



A = 



and Smoothness is defined as: 



S = 



(2) 
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where Is(i,j) is the same pixel in the image after smoothing. 

The Gini parameter is an economic indicator of equality, 
i.e., G=l; all the flux in one pixel, G=0; equal values for 
all pixels in the object. We use the implementation from 
lLotz et all l|2004 ): 
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where U is the intensity of pixel i in an flux-ordered list of 
the n pixels in the object, and I the mean pixel intensity. 
They also introduced the M20 parameter: 
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for which Y^U < 0.2 hot is true. (5) 
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where M to t = SM, = EIi[(xi -x c ) 2 + (y; -y c ) 2 ]. The center 
of the object is at (x c ,y c ). and pixel k is the pixel marking 
the top 20% point in the flux-ordered pixel-list. 

Instead of the intensity of the pixel (Ii) one can use the 
second order moment of the pixel {Mi = Ii[(xi — x c ) 2 + (yt — 
j/ c ) 2 ]) in equation^] This is our Gm parameter: 



Gm — 



M t n(n - 1) 



Hi(2i-n-l)\Mi\ 



(6) 



which is an indication of the spread of pixel values with dis- 
tance to the galaxy centre. The combination of these param- 
eters quantify the morphology of H I column density maps. 
We refer the reader for a broader discussion of these param- 
eters i n the previous papers of this series l|Holwerda et al.l 
l2011bl ich. 

As input these parameters need an estimate of the cen- 
tre of the object and a definition of the area over which they 
need to be computed. The centres of the objects are a given, 
in the case of the simulations as the centre of the dark matter 
halo. In the case of surveys as the centre of the light distribu- 
tion. To define the extent of the H I disk , we use a threshold 
of Nhi = 3 x 10 20 cm~ 2 . This defines the area over which 
the morphological parameters are computed. This thresh- 
old is the stated observat i onal limit of the W HISP survey 
jvan der Hulst et al.|[200ll ; lvan der HulstJl2002l l. and it is in 
between the two H I column density thresholds we used in 
iHolwerda et all (|2009l . l2011bh for the e xtent of the the stella r 
and gas disks in the THINGS survey (|Walter et alj|2008bj ). 

Starting with the snapshots of the H I in simulations, we 
explore in this paper the evolution of the above morpholog- 
ical parameters over time, both during a violent merger as 
well as the passive evolution in a few isolated disks. We de- 
termine how long a merging disk spends in the morphological 
parameter space that hosts mergers as a function of view- 
ing angle and physical disk characteristics. We found this 
space both from literature definitions for optical morphol- 
ogy and empirically in the WHISP sample for Hi morphol- 
ogy (|Holwerda et al.ll201ld ). Visibility time scales can then 
be used to convert an observed fraction of merging galax- 
ies into a merger rate -the number of galaxies merging in a 
given volume per Gyr- as app lied to the whole WHISP sam- 
ple in the companion paper IHolwerda et al.l (|2011al ). One 
caveat that we need to point out from the outset, is that 
our merger simulations are for 1:1 merger of massive spiral 
disks and the WHISP sample consists of a mix of galaxy 
masses (see Appendix B in the online version of the paper). 
Future suites of merger simulations should include a wider 
mass range and gas fraction. 



3 SIMULATIONS OF MERGING AND 
ISOLATED SPIRAL DISKS 

Our pri mary set of s imulat ions is a suite of equal mass merg- 
ers from lCox et "all (|2006aT ). from which an atomic hydrogen 
component was estimated by one of us (T.J. Cox) assum- 
ing thermal equilibrium. Hi maps for a control sample of 
isolated spiral disks, passively evolving for 2.5 Gyr, were 
constructed in the same manner. 

To gauge the importance of different implementations 
of a Milky Way size disk in a merger simulation, we also 
use a single merger simulation and isolated disk simulation 



Table 2. Scales of the different components of the disk galaxies. 
Component 



Cox et al 2006 Wcniger et al 2009 
(kpc) (kpc) 



Disk 

- scalclength, R^ 4.0 

- scaleheight, zq 1.0 
Gas scalelength 16.5 
Bulge scalelength, Ri, 0.45 



4.0 
0.4 



from IWeniger et all (2009), also converted to an H I column 
density map. The code, assumed physics and timescales are 
different for the Weniger simulations (see £|4.3|1 . 

All the merger simulations are for equal mass mergers 
of two large spiral galaxies. The suite from Cox et al. and the 
single simulation from Weniger et al. are for spiral galaxies 
of similar scale and mass (see Tables Q] and [2} • 

3.1 Simulations from Cox et al. 2006 



The merger simulations suite from ICox et all d2006ah use 
the N-body/SPH code G ADGET (|Springel et alj l200ll ; 



ISpringel fc Hernquistl [20021 ). It conserves entropy and fea- 
tures additional routines that track the radiative cooling of 
gas and star formation (Volker Springel, private communi - 
cation to T.J. Cox in 2001). In addition. ICox et alj (|2006aT ) 
implemented several new features themselves: stellar feed- 
back, met allicity- dependent cooling and the ability of each 
gas particle to spawn multiple new stellar particles. 

The simulations are for the isolated disks as well as the 
major mergers with variations in disk properties and view- 
ing angles (Table [3}. All these simulations run for 2.5 Gyr 
in steps of 100 Myr. Images are surface brightness values in 
10 4 h M© pc~ 2 with a pixelscale of 150 pc. The high reso- 
lution of these simulations enabled us to construct H I maps 
with sampling much finer than any current survey (See Ta- 
ble B.l [in the online appendix), and closer to the resolution 
typical for optical images. The nominal disks in these sim- 
ulations are spirals with a ~ 1O 9 M0 gas disk (unless noted 
differently), which, in the case of mergers, interacts with an 
identical disk. 

Significant disk properties that are varied are mass 
(high-m and low-m), and the treatment of the feedback from 
star-formation on the interstellar matter (ism2). We con- 
sider three orientations: face-on, at a 45° angle (orienta- 
tion2) and edge-on {edge-on). Examples of the face-on disks, 
isolated and merger, are shown in Figure [T] and Appendix A 
{online version). For more details o n th ese simulations, w e 
refer the reader to lCox et all l|2006al ) and lCox et"aH (|2006rJ ). 



3.2 The single merger simulation from Weniger et 
al (2009) 

As a consistency check, we apply o ur morphological cod e 
to the single merger simulation from IWeniger et al.l {2009), 
as well as a sing le isolated disk simulation. Following 
lHarfst et al.l (|20060 . they constructed the initial galaxies 
by firstly generating a disk/bulge/hal o - system using the 
method of iKuiiken fc Dubinskil (1995). In their next step, 
one fifth of the stellar disk particles were transformed into 
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Table 1. Mass fractions and number of particles in the different components in both simulations. 

Cox ct al 2006 Weniger et al 2009 

Component Mass fraction Nr of particles Mass fraction Nr. of particles 



Bulge-stars 

Disk-stars 

Clouds 

Gas-particles 

DM-particles 



0.012 
0.048 

0.065 
0.8 75 



10 000 
30 000 

30 000 
100 000 



0.069 
0.118 
0.016 
0.008 
0.789 



10 000 
75 098 
36 347 
9017 
100 000 



Total Mass 



8.12 xlO n M 170 000 



2.46 xlO"M 



230462 



Table 3. Cox et al Suite of Merger Simulations. 



isolated 


face-on 


reference, 3.34 X 1O 1O M0 total disk mass 


edge-on 


same but edge-on perspective 


high-mass 


face-on perspective on 8.7 X 10 10 Mq 


low-mass 


face-on perspective on 1.24 X 10 10 Mq 


mergers 


face-on 




run 1 


face-on perspective on two 




3.34 X 10 10 M Q disks collision 


run 2 


same with slightly different collision orbit. 


orientation 




45° 


same at a 45° camera angle 


edge-on 


same for edge-on camera 


disk mass 




high-mass 


face-on perspective on two 8.7 X lO lo Af0 


low- mass 


face-on perspective on two 1.24 X 1O 1O M0 


different ISM 


different ISM feedback treatment, 




see Cox et al. f2006a1 


high-mass 


with 8.7 x 1O 1O M disks. 


low- mass 


with 1.24 x 10 10 M Q disks. 



3.3 Limitations 

The suite of simu lations of isolated and merging disks form 
Cox et al ] (l2006al) as well as the implementation check from 
Weniger et alj (2009) arc limited in their scope of disk (and 
merger) properties, i.e., only a large, gas- rich spiral disk col- 
liding with its twin or evolving passively. There are only 
three different disk masses and camera angles, and no change 
in merger impact parameter of incident angle (see Table [3]). 
We limited our scope to this suite as we could easily ob- 
tain H I maps for these simulations and still be consistent in 
treatment of physics and particle resolution across the sam- 
ple (with the Weniger implementation as the exception). 
Our second reason was that simulations of isolated disks 
with identical physics and resolution evolving for a similar 
amount of time were available. 

Future simulations of the evolution of Hi morphology 
should include more minor mergers, a greater range in disk 
mass, a range in gas fraction, different impact parameters 
and incident angles of the mergers, as well as more cam- 
era angle s. The suite of simulations presented in lLotz et al.l 
l|2010al lbr) explore parameter space more, generating optical 
morphology measures. A similar effort will be needed in H I 
before volume-limited all-sky surveys in Hi are completed 
to complement the initial effort presented in this paper. 



molecular clouds and finally the diffuse ISM was added as an 
initially slowly rotating homogeneous sphere, which will col- 
lapse in the fi rst 200 Myr forming a warm disk. The collision- 
less model of lKuiiken fc Dubinskil l|l995h realises an equilib- 
rium configuration but in this case, the equilibrium is af - 
fected by adding the ISM. Therefore, I Weniger et"al] (|2009r ) 
first follow the systems evolution, until a quasi-equilibrium 
is established. The numerical integration is done by means 
of a T REE-SPH code combin ed with the sticky particle 
method (|Theis fc He nslcr 1993|). G ravitational fo rces are de- 
termined by the DEHNEN-Tree l|DehnerJ I2OO2T ) . The cold 
gas (T < lO 4 -^) in this simulation was converted to the H I 
map, assuming that the fraction of the cold gas which is 
observable as Hi, depends on the thermal pressure. The HI 
map was smoothed by a gaussian with a FWHM of 1 kpc. 
The centre of the interacting disk is taken to be the centre of 
the bulge stars distribution. Time steps are 10 and 50 Myr 
for the isolated disk and merger sim ulation, respectively . For 
more details we refer the reader to lHarfst et all (|2006u and 
IWeniger et all (|2009r ). 



3.4 Shaping the Hi Morpholgy 

The morphology of the H I disk is influenced by many differ- 
ent processes, only one of which is gravitational interaction 
with a neighbouring galaxy. The morphology is also deter- 
mined by the level of star-formation (triggered or not by the 
interaction) , which consumes gas, and the stellar winds and 
supernovae feedback that create the typical "effervescent" 
look of Hi maps. An additional effect from star-formation 
(or an active AGN) can be that much of the gas in a galaxy 
is heated or ionized, and thus no longer visible in the H I 
map. Depending on the temperature and dust fraction of 
the ISM, much of the higher density parts of the spiral gas 
disk will be in the form of molecular hydrogen, also remov- 
ing it from the Hi map. Ram-pressure effects experience 
during travel through the intergalactic medium of a cluster 
may also contribute significantly to the appearance of the 
Hi disk. 

The simulations treat ISM feedback from star-formation 
comprehensively but for now, we set aside the effects of 
an AGN switching on or the interaction happening in a 
dense cluster medium . In a further paper in this series 
l|Holwerda et al.ll2011dh . we do explore the how much ram- 
pressure and H I stripping of galaxies change these parame- 
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ters usi ng the VLA Imagi ng of Virgo spirals in Atomic gas 
fVIVA. IChung et alJlioog ). We find that stripped Hi disks 
are easily identified with low values of Concentration. The 
effect of moderate ram-pressure on a sample of Hi disks 
will, however, be more difficult to quantify. The phenomena 
of Hi lopsidedness, which may be related to ram-pressure 
stripping see ms to have only a mode rate effect on these pa- 
rameters (see iHolwerda et al.ll2011ch . 

Therefore, we compare morphological evolution of H I 
disks with the same physics (feedback, ISM heating, etc.), 
in identical disks, and only change one critical aspect: in 
some simulations these disks merge with a similar disks, in 
others, they evolve passively. 



4 H i MORPHOLOGY OF MAJOR MERGERS 

We have run four isolated d isks and eight m erger simula- 
tions over a time of 2.5 Gyr (|Cox et al.ll2006ah and a single 
simulation of a maj or merger and one is olated face-on disk 
that ran for 3 Gyr (|Weniger et al.ll2009l ). The resulting ta- 
bles of morphological parameters over time can be found in 
Appendix B (online version) for both isolated and merging 
disks. The default perspective is face-on, with variation in 
disk mass and treatment of the ISM (Table [3j. 



4.1 Merging and Isolated Disks 

Figure [2]shows the progression over time of three parameters 
for both the merger scenarios (black lines) as well as the 
isolated disks (gray lines) for the Hi disks. We note how 
these three parameters already separate isolated disks and 
the merging ones as well as identify the times the galaxies 
are interacting the strongest. 

All three parameters show a clear parameter space 
where the galaxy disk is almost certainly merging (Gm > 
0.6, M20 ~ —3, and A > 0.5) and a part where it is unlikely 
to be merging. In the case of M20 and Gm, the merging disks 
return to the isolated disk's value periodically. Their values 
spike at the time of close passage of the merging companion. 
Variation of the disk mass (low-m and high-m simulations) 
do not change the overall morphological signature apprecia- 
bly. 



4.2 Edge-on Disks 

Edge-on galaxies are morphologically the least disturbed 
perspective on a merger. Figure [3] shows the isolated disk 
and merging disk in cold gas, just from the edge-on per- 
spective (dotted lines in Figure [5}. The merger is still mor- 
phologically selected by our selection criterion (see £ 14.41) . 
albeit not very often. However, in the case of the edge-on 
perspective, any signature in the dynamical profile of the 
galaxy will be the strongest, i.e., there will be a clear devi- 
ation of the typical "double-horned" profile of the Hi line. 
Thus, our morphological approach and a kinematic approach 
to selecting interacting spirals in large Hi survey are very 
complementary. 



4.3 Comparison to Weniger et al 

Figure [4] shows the values of Asymmetry, M20 and Gm of 
the H 1 maps for a merging and an iso lated galaxy from the 
simulation from lWeniger et al.l (|2009l ). Asymmetry shows a 
steady rise with time for the isolated disk and clear spikes 
during the major encounters. M20 displays a very similar 
evolution as the plot in Figure[2]with some notable evolution 
by the isolated disk. Gm spikes at the close encounters, very 
similar to Asymmetry. 

The morphological selection criteria (see next section), 
perform well on the Weniger et al. simulation, selecting the 
merging galaxy for part of the time and mostly not selecting 
the isolated galaxy's Hi disk. The times, in billions of years, 
the Weniger et al merger are selected by the various mor- 
phological criteria are listed at the bottom of Table U The 
general similar behavior of the Weniger et al. simulation to 
the suite from Cox et al. is an indication of the robustness 
of the morphological selection. 

4.4 Merger Selection Criteria 

The goal of this paper is to obtain a typical time scale that 
mergers are visible as disturbed Hi morphologies. The visi- 
bility time depends on the selection criteria, observed wave- 
length and orientation of the merger. We can now explore 
how the visibility times for each visibility criterion relate to 
the different physical properties of each merger run, e.g., gas 
mass, viewing angle, etc. 

Originally these parameters were envisaged to classify 
the morphologies of galaxies and hence a single parameter 
criterion will not necessarily work. Interacting spiral disks 
are a subset of t he pa rameter-space occ upied by late-types. 
IConselicd (120031 ) and lLotz et all (|2004f ) introduced several 
different criteria for the selection of merging syste ms in their 
respec tive parameter systems. For optical data, IConselicel 
(2003) define the following criterion: 

A > 0.38, (7) 

with some authors requiring A > S as well. In general, highly 
asymmetric galaxies are candidate merging systems (Figure 

lLotz et al. I l|2004h added two different criteria using Gini 
and M20: 

G > —0.115 x A/20 + 0.384 (8) 
and 

G > -0.4 x A + 0.66 or A > 0.4. (9) 

The latter being a refinement of the Conselice et al criterion 
in equation [7] 

These criteria were developed for optical morphologies, 
typically observed Johnson-B or sdss-g. Typical optical spa- 
tial resolution (~ 2") is a factor of a few higher than the 
typical spatial resolution of the Hi maps (6-12", depending 
on the survey). Theref ore, in the third paper in this series 
l|Holwerda et al.ll201ld) . we defined several possible criteria 
specifically for the Hi perspective using the CAS-G/A/20- 
Gm space of the WHISP survey H imap sample. We defined 
the Gini parameter of the second order moment, Gm and a 
criterion that selected most interacting galaxies: 

G M > 0.6, (10) 
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Figure 1. Three Hi snapshots (at times T=0. 1,1.2 and 2.4 Gyr) of two of our simulations, the nominal isolated, passively evolving spiral 
disk (left) and the interacting and one instance of the merging disk (right). Contours are at 0.5, 1.0, and 2.5 Mqpc~ 2 to highlight the 
extent of the disks at low column densities. 
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T (Gyr) T (Gyr) T (Gyr) 

Figure 2. Asymmetry, M20 and Gm as a function of time in the isolated (gray) and mergers (black) for the Hi disks from lCox et al.l 
(200630). The disks are differentiated between nominal disks (solid lines), higher mass (short dashed lines), lower mass (long dashed 
lines), with different ISM treatment (high- mass; dot-short dash, low- mass; dot-long dash), and edge-on (dotted lines). The thick lines 
are the times A < —0.2 X M20 + 0.25, one of our selection criteria. Isolated disks (gray lines) have less Asymmetry (A < 0.5), high values 
for M20 (close to 0) and a low Gm values. Based on these timelines, the mergers do occasionally return to quiescent parameter space 
but differentiate well for a large fraction of the time. 




T (Gyr) T (Gyr) T (Gyr) 

Figure 3. Edge-on: Asymmetry, M20 and Gm of the edge-on cases as a function of time in the isolated (gray) and mergers (black) of 
the cold gas. The thick lines are the times A < —0.2 X M20 + 0.25, one of our selection criteria. 
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Figure 4. Asymmetry, M20 and Gm of both galaxies in the simulation of Wcniger ct al. The morphology of the merging galaxy 

in H I (black lines) and isolated galaxy (gray lines) are shown as a function of time. As an example, the thick black lines are those times 
of the simulation the merger would be selected by the A < —0.2 X M20 + 0.25 criterion. 
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a criterion that seems to be corroborrated by Figure [2] 

Earlier in this series, we speculated that a combina- 
tion of Asymmetry and M20 co uld well be used to se lect 
interaction in H I morp hology in |Holwerda et al.ll2011bl ) . In 
iHolwerda et~ai] |201ld) . we defined this criterion as: 

A > —0.2 x M20 + 0.25. (II) 

Finally, we also defined one b ased on Concentra tion and 
M20, following the example of the lLotz et al. I Hil) criteria 
(eq.|5]and|!|l): 

C > -5 x M 20 + 3. (12) 

In l|Holwerda et al.ll2011bl ), we found that this last criterion 
selected to both the correct fraction of interacting galaxies 
as well as agree most often with the previous visual identi- 
fications in the case of individual WHISP galaxies. 

We can now explore the timescales of each of these cri- 
teria, both those from the literature as well as those we 
determined in the WHISP sample subset. The most reli- 
able merger selection criterion would be the one that se- 
lects the right fraction of galaxies out of a given sample 
for a long timescale, with little contamination from isolated 
galaxies. Preferably, viewing angle or physical disk charac- 
teristics (e.g., gas mass) do not influence the timescale overly 
much. A long selection time scale would subsequently ensure 
the more accurate estimate of the merger rate from a given 
volume of galaxies. 

Figure [5] shows the relations between the morphological 
parameters for the Hi gas disks in the Cox et a l. simulations, 
similar to Figure 4 in IHolwerda et al.1 l|2011ch . from which 
we determined the selection criteria for H I morphologies (eq. 
I10fll2|l . The literature and our merger selection criteria are 
marked with dashed and dotted lines respectively. The times 
an Hi disk is selected by the various selection criteria are 
listed in Table H 

Figure[6]shows the six criteria feq.l7l fT2]) as a function of 
time for all our simulations, both isolated and merging disks. 
If the y-axis values are positive, the merger would have been 
selected by this criterion. Immedia tely, it is apparent, that 
the three criteria that we defined in IHolwerda et all (|201ld ) 
perform very well, separating merging galaxies from isolated 
disks. It also is clear that selection happens most often in 
the fist 1.5 Gyr of the merger simulations and that some of 
the criteria could be adjusted for the native resolution of 
these simulations. 

The selection by Asymmetry alone (eq. [7| or similar) 
does not translate well to Hi. Asymmetry alone selects iso- 
lated and merging galaxy almost equal times (Figure [6] 
top left panel, Table |4| . This validates what we found for 
this criterion in the W HISP sample and our suspicion in 
IHolwerda et~aH (|2011bh that H 1 Asymmetry is influenced by 
other factors than merging. From Figure [6] one could con- 
clude that a more strickt Asymmetry criterion would work 
better but this is only valid for the resolution of our sim- 
ulations. The Gini-M2o criterion (eq. [8} would need to be 
modified for the H I perspective and performs poorly in the 
separation of mergers from isolated disks in its current def- 
inition (Figure [6] top middle panel) . The Gini- Asymmetry 
criterion (eq. [9] Figure [6] top right panel) performs much 
better, very rarely selecting the isolated disk and selecting 
mergers on average for 0.44 Gyr out of the 2.5 Gyr the sim- 
ulations ran. 



However, these simulations ran on the maximum spatial 
resolution (150 pc sampling). The selection criteria perform 
a little poorer for lower resolution (Appendix B in the elec- 
tronic edition). For example, the G-A selection criterion (eq. 
[SJ, still performs well for the WHISP survey observational 
parameters (Table ??) but with a typical selection time scale 
of 0.15 Gyr. These merger visibility times are close to the 
ones typically quoted for these sele ction criteria app lied to 
optical images of stellar disks fe.g.. lLotz et al.ll2010al lbh. 

Of the literature criteria, only the G-A one appears to 
translate directly from optical to H I morphology. Our selec- 
tion criteria (eq. I10H12[) were based on a WHISP subsample 
of H 1 column density maps for which we had visual estimates 
of interaction. Our Gm criterion (eo. HOp cleanly select merg- 
ers for 0.69 Gyr of the 2.5 Gyr runtime (Figure(6] bottom left 
panel). The selection becomes diluted, however, when the 
simulations are smoothed to the WHISP resolution (Table 
B.3 online appendix) or worse (e.g., VIVA, Table B.4 online 
appendix). The A-M20 selection criterion (eq. Ill|l appeared 
to work extremely well in the WHISP data. Yet, in the simu- 
lations, it selects isolated disks more often, on average, than 
merging ones (Figure[6] bottom middle panel). One can con- 
ceivable still use criteria like this one, provided the selected 
galaxy fraction is corrected for the contamination by iso- 
lated disks. Or the criterion is adjusted to the resolution of 
the data. The Concentration- M20 criterion (eq . 112(1 performs 
better, selecting mergers for 1.19 Gyr versus isolated disks 
for 0.69 Gyr, but would also need a substantial correction 
for the isolated disk contaminations (Figure [6] bottom right 
panel). We note that the main reason isolated disks are se- 
lected by this criterion are the low gas-mass and edge-on disk 
simulation and this criterion performs much better for the 
low inclination disks. Contamination increases again when 
the simulations are degraded to WHISP resolution (Table 
B.3 in the online appendix). 

In IHolwerda et~ai1 (|201ld ). we found that the 
Concentration- M20 criterion flagged not only the correct 
fraction of galaxies as merging but also agreed in many cases 
with the visual classifications on which galaxies were inter- 
acting. Here we find that the merger visibility time scale is 
in a similar range or even a little better than the selection 
criteria from the literature for optical images of disks (typ- 
ically ^ 1 Gyr). Thus, depending on the volume surveyed 
and the spatial resolution of the survey, an Hi survey of 
disks can provide an estimate of the merger rate at least as 
accurate as any optical one using these parameters. 

lLotz et all l|2010af) note that gas-rich mergers tend to be 
selected more often in the optical morphological selections, 
because gas rich mergers trigger more star-formation, which 
leaves brighter tidal features in the disks before full merger. 
A similar selection biases is pertinent for Hi morphology: 
more gas-rich disks will stand out in 21 cm. observations. 
Lowering or increasing the gas-mass of the disks, does influ- 
ence the selection times, as does a different treatment of the 
ISM physics (ism2 or the Weniger simulation), increasing 
the variance in merger visibility times. 

The perspective or camera angle on the merger does in- 
fluence the selection time for all the different criteria. Oddly, 
the edge-on perspective appears to occasionally raise the 
time a merger gets selected by some of the morphological 
parameter criteria. However, this can be expected to be the 
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poorest viewing angle, regardless whether or not the mor- 
phology is determined in Hi or optical. 

The changes in perspective, gas disk mass and ISM 
treatment for different simulations give an indication of the 
spread in merger visibility times in Table 2] which is sub- 
stantial. Therefore, we can determine an accurate merger 
fraction for a sample of galaxies associated with a known 
volume, yet the merger rate will still be subject to some un- 
certainty due to the range in possible merger visibility times. 
With these average values in Table [4] in hand, morphologi- 
cally selected fraction of an H I survey can n ow be converted 
to a m erger rate with reasonable accuracy l|Holwerda et al.l 
l2011ah . 



5 CONCLUSIONS 

In this paper we obtained an understanding of the morphol- 
ogy of Hi disks during a major merger. Observationally, we 
have already found that the morphological parameters, cur- 
rently in use to classify galaxies and identify mergers in op- 
tical and UV wavelengths can be applied very productively 
to Hi maps. Applying these to the cold gas and Hi maps 
based on N-body simulations, we conclude the following: 

1. The current merger simulations do an excellent job re- 
producing the stellar disks and its morphology as well as 
an accurate cold atomic hydrogen (Hi) map. Morphology of 
the H I maps follows those seen in observations well both in 
appearance (Figure[TJ and quantified with CAS, G and M20 
(Figure [5j). 

2. Asymmetry, M20 and Gm are good parameters to dis- 
tinguish mergers from isolated disks (Figure [2] and [4} but a 
combination of morphological parameters often works better 
(Figure El). 

3. Both the Weniger et al and the Cox et al. merger simu- 
lations show very similar behavior in morphology (Figure [2] 
and [4}. Different implementations of the ISM physics need 
not result in dramatically different global H I morphology. 

4. Disk gas mass, orientation or the treatment of ISM 
physics do not substantially change the separation in pa- 
rameter space of merging and isolated disks in these three 
parameters (Figure |2J. 

3. Even the edge-on perspective on a merger produces 
a morphological signature in the Hi column density maps 
(Figure[3]), yet this remains the poorest viewing angle, often 
confusing merging and isolated disks. 

4 Of the merger criteria from the literature, the one based 
on the Gini and Asymmetry performs well for the Hi per- 
spective, selecting mergers for 0.4 Gyr out of a runtime of 
2.5 Gyr. 

5 The merger criteria defined in iHolwerda et all (|201ld ) 
specifically for Hi show a range of merger selection time 
scales as well as some contamination from isolated disks 
(Figure [S] and Table 3}. Of these three, the Gm selection 
appears to be the cleanest selection with 0.69 Gyr visibility 
time. 

6 The criteria we found to perform the best on the 
WHISP sample, C > —bA'ho + 3, has a visibility time of 
~ 1 Gyr. but with some contamination by isolated disks, 
mostly low gas-mass and edge-on ones. 



In the future, more detailed and comprehensive suites 
of simulated H I disks will become available, which will allow 
for determinations of merger time scales for both dynami- 
cal and morphological signatures in H I observations. Such 
simulations will become necessary when the large volume 
surveys on the SKA precursors commence. Notably, the all- 
sky Hi survey, that will result from the WALLABY project 
with the ASKAP radio telescope and the Northern H I sur- 
vey with the APERTIF instrument on WSRT combined, will 
produce hundreds of thousands of Hi data cubes, at similar 
resolutions as WHISP. Automated merger classification of 
this surveys will be the only feasible way to delineate sam- 
ples, as well as determine a more accurate merger rate in the 
local Universe from an H I perspective. Hopefully, this obser- 
vational effort will be accompanied by an expanded suite of 
merger simulations to obtain a more accurate mean selection 
rate and possibly a spread due to difference in IS M, orienta- 
tion an d disk masses. Following the reasoning of lLotz et al.l 
l|2010al lbh. we will need to determine merger visibility time 
scales for the gas-rich and unequal mergers for which these 
upcoming surveys will be especially sensitive. 

Another possible future use of H I morphological param- 
eters is to combine them with at other wavelengths, e.g., 
optical. The different times at which the disk becomes dis- 
turbed -earlier in Hi, followed by optical- could be used to 
estimate in what phase any given merger is in. 
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1. A > 0.4 from lConseiic3 l|2003h 



2. G > -0.115M20 + 0.38 4 fromlLotz et alJ J2004h 

3. G > -OA A + 0. 66 fromlLotz et all \2004) 

4. G M > 0.6 from lHolwerda et alJ <2011ch. 

5. A < -0.2M 20 + 0.25 fromlHolwerda et alJ j20lT3 'l. 

6. C> -5M 20 + 3 from lHolwerda et al.l j201ljT 

7. The combination of G M > 0.6 and A < -0.2M 20 + 0.25. 

8. The combination of A < -0.2Af 20 + 0.25 and C > -5M 20 + 3. 

Table 4. The times in Gyr an H I simulation from lCox et al. ( 2006a) is selected by the different criteria. Criteria arc from Holwerda ct al. 
and those adopted from the literature. 
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